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We characterize the complexity of geodesic paths on a curved statistical manifold A4 S through the 
asymptotic computation of the information geometric complexity Vm s an d the Jacobi vector field 
intensity Jjvi s - The manifold M 3 is a 2Z-dimensional Gaussian model reproduced by an appropriate 
embedding in a larger 4Z-dimensional Gaussian manifold and endowed with a Fisher-Rao information 
metric (O) with non-trivial off diagonal terms. These terms emerge due to the presence of a 
correlational structure (embedding constraints) among the statistical variables on the larger manifold 
and are characterized by macroscopic correlational coefficients ru- First, we observe a power law 
decay of the information geometric complexity at a rate determined by the coefficients and 
conclude that the non-trivial off diagonal terms lead to the emergence of an asymptotic information 
geometric compression of the explored macrostates on M a . Finally, we observe that the presence 
, of such embedding constraints leads to an attenuation of the asymptotic exponential divergence of 

Q ■ the Jacobi vector field intensity. 

PACS numbers: Probability Theory (02.50.Cw), Riemannian Geometry (02.40.Ky), Chaos (05. 45. -a), Com- 
l/-) . plexity (89.70.Eg), Entropy (89.70.Cf). 

S 1 : 

O-t- I. INTRODUCTION 

Characterizing and understanding the mystery of the origin of life and the unfolding of its evolution are perhaps the 
leading arguments motivating the quantification of the extremely elusive concept of complexity [U-IH . Of course there 
i " i ' are more pragmatic reasons that justify the study of complexity, for example the problem of quantifying how complex 
is quantum motion. This issue is of primary importance in quantum information science. However, our knowledge of 
the relations between complexity, dynamical stability, and chaoticity in a fully quantum domain is still not satisfactory 
[1, The concept of complexity is very difficult to define, its origin is not fully understood [6rJTT] and it is mainly 
for this reason that several quantitative measures of complexity have appeared in the scientific literature GLU- In 
classical physics, measures of complexity are settled in a better way. The Kolmogorov-Sinai metric entropy [12 1, that 
ly-^ . is the sum of all positive Lyapunov exponents [l3j , is a powerful indicator of unpredictability in classical systems and 
it measures the algorithmic complexity of classical trajectories (l3 - [l7| . Other known measures of complexity are the 
logical depth [l^, the thermodynamic depth the computational complexity [2(|, the stochastic complexity [2l[ 
and many more. Ideally, a good definition of complexity should be mathematically rigorous and intuitive at the same 
time so that it allows to tackle complexity-related problems in computation theory and statistical physics as well. Of 
course, a quantitative measure of complexity is truly useful if its range of applicability is not limited to few unrealistic 
applications. It is also for this reason that in order to properly define complexity measures, it should be clearly stated 
the reasons why one is defining such a measure and what it is intended to capture. 

It is known that classical complex systems exhibit local exponential instability and are characterized by positive 
Lyapunov exponents [22j . Furthermore, the study of the relationship between entropy and the complexity of trajec- 
tories of a dynamical system has always been an active field of research [ll-QJ]. In particular, in [14J it was shown 
that the algorithmic complexity of trajectories of points in a dynamical system is asymptotically equal to the entropy 
of the system. 

It is commonly accepted that one of the major goals of physics is modeling and predicting natural phenomena using 
relevant information about the system of interest. Taking this statement seriously, it is reasonable to expect that the 
laws of physics should reflect the methods for manipulating information. Indeed, the less controversial opposite point 
of view may be considered where the laws of physics are used to manipulate information. This is exactly the point of 
view adopted in quantum information science where information is manipulated using the laws of quantum mechanics 

[H. 

Here we make use of the so-called Entropic Dynamics (ED) [24] and Information Geometrodynamical Approach 
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to Chaos (IGAC) [25|, |26[. ED is a theoretical framework that arises from the combination of inductive inference 
(Maximum Entropy Methods, [27]-[3(|) and Information Geometry [3l|. The most intriguing question being pursued 
in ED stems from the possibility of deriving dynamics from purely entropic arguments. This is clearly valuable 
in circumstances where microscopic dynamics may be too far removed from the phenomena of interest, such as 
in complex biological or ecological systems, or where it may just be unknown or perhaps even nonexistent, as in 
economics. The applicability of ED has been extended to temporally-complex (chaotic) dynamical systems on curved 
statistical manifolds and relevant measures of chaoticity of such an information geometrodynamical approach to chaos 
have been identified [25[. IGAC arises as a theoretical framework to study chaos in informational geodesic flows 
describing physical, biological or chemical systems. A geodesic on a curved statistical manifold Ais represents the 
maximum probability path a complex dynamical system explores in its evolution between initial and final macrostatcs. 
Each point of the geodesic is parametrized by the macroscopic dynamical variables {0} defining the macrostate of the 
system. Furthermore, each macrostate is in a one-to-one correspondence with the probability distribution {p(A|0)} 
representing the maximally probable description of the system being considered. The quantity X is a microstate of 
the microspace X . The set of macrostates forms the parameter space T>q while the set of probability distributions 
forms the statistical manifold Ais- IGAC is the information geometric analogue of conventional geometrodynamical 
approaches |32l |33| where the classical configuration space Te is being replaced by a statistical manifold Ais with 
the additional possibility of considering chaotic dynamics arising from non conformally flat metrics (the Jacobi metric 
is always conformally flat, instead). It is an information geometric extension of the Jacobi geometrodynamics (the 
geometrization of a Hamiltonian system by transforming it to a geodesic flow [34|). The reformulation of dynamics 
in terms of a geodesic problem allows the application of a wide range of well-known geometrical techniques in the 
investigation of the solution space and properties of the equation of motion. The power of the Jacobi reformulation 
is that all of the dynamical information is collected into a single geometric object in which all the available manifest 
symmetries are retained- the manifold on which geodesic flow is induced. For example, integrability of the system 
is connected with existence of Killing vectors and tensors on this manifold. The sensitive dependence of trajectories 
on initial conditions, which is a key ingredient of chaos, can be investigated from the equation of geodesic deviation. 
In the Riemannian [32j and Finslerian [33[ (a Finsler metric is obtained from a Riemannian metric by relaxing the 
requirement that the metric be quadratic on each tangent space) geometrodynamical approach to chaos in classical 
Hamiltonian systems, an active field of research concerns the possibility of finding a rigorous relation among the 
sectional curvature, the Lyapunov exponents, and the Kolmogorov-Sinai dynamical entropy [35j . 

In this article, inspired by the work presented in (36j . using statistical inference and information geometric tech- 
niques, we characterize the complexity of geodesic paths on a curved statistical manifold M s through the asymptotic 
computation of the information geometric complexity Vm, and the Jacobi vector field Jm s - The manifold A4 S is a 
2/-dimcnsional Gaussian model reproduced by an appropriate embedding in a larger 4^-dimensional Gaussian manifold 
and endowed with a Fisher-Rao information metric g M „ (O) with non-trivial off diagonal terms. These terms in the 
information metric on the embedded manifold emerge due to the presence of a correlational structure (embedding 
constraints) among the statistical variables on the larger manifold and are characterized by macroscopic correlational 
coefficients r^. First, we observe a power law decay of the information geometric complexity Vm s ( r ) ~ exp [Sm 3 i T )] 
at a rate determined by the coefficients and conclude that the non-trivial off diagonal terms lead to the emergence 
of an asymptotic information geometric compression of the explored macrostates on Ai s - Finally, we observe that 
the presence of such embedding constraints leads to an attenuation of the asymptotic exponential divergence of the 
Jacobi vector field intensity. 

The layout of this article is as follows. In Section II, we present few remarks on the theoretical structure of the 
IGAC and outline few selected applications concerning the complexity characterization of geodesic paths on curved 
statistical manifolds. In Section III, we describe the 2^-dimensional curved statistical model considered, the embedded 
Gaussian model endowed with a Fisher-Rao information metric with non-trivial off diagonal terms. In Section IV, we 
present the asymptotic computation of the information geometric entropy <Sm s ( t )- We observe a power law decay 
of the information geometric complexity Vm b ( t ) at a rate determined by the coefficients and conclude that non- 
trivial off diagonal terms lead to the emergence of an asymptotic information geometric compression of the explored 
macrostates on the statistical configuration manifold considered. In Section V, we present the asymptotic computation 
of the Jacobi fields on X5. We observe that the presence of the embedding constraints lead to an attenuation of the 
asymptotic exponential divergence of the Jacobi field intensity. Finally, in Section VI we present final remarks and 
suggest further research directions. 

II. ON THE IGAC: REMARKS AND APPLICATIONS 

In this Section, we present few remarks on the theoretical structure of the IGAC and outline few selected applications 
concerning the complexity characterization of geodesic paths on curved statistical manifolds. A more detailed review 
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appears in [25J. 

A. Remarks 

As stated in the Introduction, the IGAC arises as an information geometric framework to study chaos and complexity 
in informational geodesic flows describing physical systems. A geodesic on a curved statistical manifold represents the 
maximum probability path a complex dynamical system explores in its evolution between initial and final macrostates 
Qi and 9/, respectively. Each point of the geodesic on a 4Z-dimensional statistical manifold represents a macrostate 
parametrized by the macroscopic dynamical variables = 04i) defining the macrostate of the system. 

Furthermore, each macrostate is in a one-to-one correspondence with the probability distribution P (X\Q) representing 
the maximally probable description of the system being considered. The set of macrostates forms the parameter space 
while the set of probability distributions form the statistical manifold. In what follows, we schematically outline the 
main features underlying the construction of an arbitrary form of entropic dynamics. First, the microstates of 
the system under investigation must be defined. For the sake of reasoning, we assume the system is characterized 
by an 2/-dimensional microspace with microstates X = (x\,..., X2i)- The main goal of an ED model is that of 
inferring " macroscopic predictions''' in the absence of detailed knowledge of the microscopic nature of the arbitrary 
complex systems being considered. More explicitly, by "macroscopic prediction" we mean knowledge of the statistical 
parameters (expectation values) of the probability distribution function that best reflects what is known about the 
system. This is an important conceptual point. The probability distribution reflects the system in general, not the 
microstates. Once the microstates have been defined, we then select the relevant information about the system. In 
other words, we have to select the macrospace of the system. 

In general, one is given a manifold of probability distributions arising from the maximum entropy formalism where 
distributions arise from the maximization of the logarithmic relative entropy subjected to some statements concerning 
averages [information constraints). Given the manifold of probability distributions, the {direct) problem is to find 
the corresponding Fisher-Rao information metric. However, not all probability distributions are generated in this 
way. For instance, probability distributions may emerge as a result of a change of variable technique (parametric 
transformation law) [37| . Furthermore, in order to do physics, we are also concerned with the following [inverse) 
problem: we want to design statistical manifolds with appropriate geometries [38j . 

B. Applications 

In the following, we outline few selected applications concerning the complexity characterization of geodesic paths 
on curved statistical manifolds. 

1. Gaussian Statistical Models in the Absence of Correlations 

In (36[, we apply the IGAC to study the dynamics of a system with I degrees of freedom, each one described by 
two pieces of relevant information, its mean expected value and its variance (Gaussian statistical macrostates). This 
leads to consider a statistical model on a non-maximally symmetric 2/-dimensional statistical manifold M s . It is 
shown that M. s possesses a constant negative scalar curvature proportional to the number of degrees of freedom of 
the system, TZm 3 = ~l- It is found that the system explores statistical volume elements on M s at an exponential rate. 
The information geometric entropy <Sm s increases linearly in time (statistical evolution parameter) and, moreover, 
is proportional to the number of degrees of freedom of the system, <Sm s r ^°° IXt where A is the maximum positive 
Lyapunov exponent characterizing the model. The asymptotic linear information geometric entropy growth may be 
considered an information-geometric analogue of the von Neumann entropy growth introduced by Zurek-Paz, a quan- 
tum feature of chaos. The geodesies on M. s are hyperbolic trajectories. Using the Jacobi-Levi-Civita (JLC) equation 
for geodesic spread, we show that the Jacobi vector field intensity Jm s diverges exponentially and is proportional 
to the number of degrees of freedom of the system, Jm 3 r ~'°° Zexp(Ar). The exponential divergence of the Jacobi 
vector field intensity Jm s is a classical feature of chaos. Therefore, we conclude that 1Zm s = ~U Jm s t ^ j ° c /exp (At) 
and Sm s t ~'°° IXt. Thus, 1Zm s , Sm s an d Jm b behave as proper indicators of chaoticity and are proportional to the 
number of Gaussian-distributed microstates of the system. This proportionality, even though proven in a very special 
case, leads to conclude there may be a substantial link among these information geometric indicators of chaoticity. 
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2. Gaussian Statistical Models in the Presence of Correlations between Microvariables 

In [39l j . we apply the IGAC to study the information constrained dynamics of a system with I — 2 microscopic 
degrees of freedom. As working hypothesis, we assume that such degrees of freedom are represented by two correlated 
Gaussian-distributed microvariables characterized by the same variance. We show that the presence of microcor- 
relations lead to the emergence of an asymptotic information geometric compression of the statistical macrostates 
explored by the system at a faster rate than that observed in absence of microcorrelations. This result constitutes 
an important and explicit connection between (micro)-correlations and (macro)-complexity in statistical dynamical 
systems. The relevance of our finding is twofold: first, it provides a neat description of the effect of information 
encoded in microscopic variables on experimentally observable quantities defined in terms of dynamical macroscopic 
variables; second, it clearly shows the change in behavior of the macroscopic complexity of a statistical model caused 
by the existence of correlations at the underlying microscopic level. 

3. Ensemble of Random Frequency Macroscopic Inverted Harmonic Oscillators 

In [38|, we explore the possibility of using well established principles of inference to derive Newtonian dynamics 
from relevant prior information codified into an appropriate statistical manifold. The basic assumption is that there 
is an irreducible uncertainty in the location of particles so that the state of a particle is defined by a probability 
distribution. The corresponding configuration space is a statistical manifold the geometry of which is defined by the 
Fisher-Rao information metric. The trajectory follows from a principle of inference, the method of Maximum Entropy. 
There is no need for additional "physical" postulates such as an action principle or equation of motion, nor for the 
concept of mass, momentum and of phase space, not even the notion of time. The resulting "entropic" dynamics 
reproduces Newton's mechanics for any number of particles interacting among themselves and with external fields. 
Both the mass of the particles and their interactions are explained as a consequence of the underlying statistical 
manifold. 

Following this line of reasoning, in [HI, 53] we present an information geometric analogue of the Zurek-Paz quantum 
chaos criterion in the classical reversible limit. This analogy is illustrated by applying the IGAC to a set of n- uncoupled 
three-dimensional anisotropic inverted harmonic oscillators characterized by a Ohmic distributed frequency spectrum. 

4- IGAC of regular and chaotic quantum spin chains 

In [ill . |42| , we study the entropic dynamics on curved statistical manifolds induced by classical probability distribu- 
tions of common use in the study of regular and chaotic quantum energy level statistics. Specifically, we propose an 
information geometric characterization of chaotic (integrable) energy level statistics of a quantum antiferromagnetic 
Ising spin chain in a tilted (transverse) external magnetic field. We consider the IGAC of a Poisson distribution 
coupled to an Exponential bath (spin chain in a transverse magnetic field, regular case) and that of a Wigner-Dyson 
distribution coupled to a Gaussian bath (spin chain in a tilted magnetic field, chaotic case). Remarkably, we show 
that in the former case the IGE exhibits asymptotic logarithmic growth while in the latter case the IGE exhibits 
asymptotic linear growth. 

III. THE MODEL 

In this Section, we emphasize the main reasoning underlying the choice of the new proposed statistical model and 
study its information geometric properties. 

A. Motivations 

We want to make reliable macroscopic predictions when only partial knowledge on the micro-structure of a system 
is available. As stated in the Introduction, the complexity of such predictions is quantified in terms of the IGE and 
the Jacobi field intensity. In this manuscript we seek an answer to the following question: isn't simpler to make 
macroscopic predictions when the available pieces of information are not independent? Stated otherwise, does an 
increase in the correlational structure of the dynamical equations for the statistical variables labelling a macrostate of 
a system imply a reduction in the complexity of the geodesic paths? It is reasonable to expect that the emergence of 
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a correlational structure in the form of constraints among the variables labelling the macrostates of a system would 
lead to a highly constrained dynamics and, consequently, to a reduction in the complexity of making macroscopic 
predictions. In what follows, we attempt to give a quantitative answer to the above-mentioned questions. 

We propose to compare the complexity of making predictions in two different scenarios. In the first scenario (Al- 
dimensional larger Gaussian model) , we assume a system with 21 degrees of freedom x k , each one being Gaussian 
distributed. The probability distribution describing the whole system is given by, 



P(X\Q) = l[p(x k \ f i k , a k ), 



(1) 



fe=i 



where X = (x 1 ,..., x 2 i), = (Md-j ^2h °"i>-j &2i) andp(x fe |^ fe , er fe ) is defined as, 



P{XkW, <Jk) 



dcf 



y/27 



■. exp 



2*2 



(2) 



In [36|, |44[ , we provided an analytical estimate for the complexity of geodesic paths on the statistical manifolds of 
distributions given in ([1}. In the second scenario (2/-dimensional embedded Gaussian model), we consider the very 
same system with 21 degrees of freedom x k - However, the microvariables x k are described by probability distributions 
characterized by statistical variables subject to a set of 2/-constraints, 

<?2j = °~2j-i and, fi 2 j = ^2j {^2j-iJ a 2j-i) (embedding constraints), (3) 

with j = I,..., I. Therefore, the probability distribution describing the whole system in this second scenario becomes, 

P (X\G) = Y[p(x 2 j-u x 2j \n 2j _ 1 , <7 2 j-i) , (4) 
where X = (x x ,..., x 2i ), = {n x , ft,-, <ri, era,—, 0"2I— l) and p (x 2 j-i, x 2j \^ 2j _ 1 , <J 2 j-i) is defined as, 



p(x 2 j-l, af2j|/i 2 j-li < T 2j-l) 



dcf 



exp 



(X2j-1 



M 2 



) 2 + [x2j - n 2 j (fhj-i, °ji-i)T 



(5) 



with j — 1,..., Our purpose is computing the complexity of geodesic paths on the 2Z-dimensional statistical 



manifold M. 



(embedded) 



= Ms = 



{P(x\e)} 



and compare it with the one obtained on the 4Z-dimensional manifold 

^Oarger) = {p( X \Q)}. We expect that the emergence of the correlational structure defined in (J3J between pairs of 
macroscopic statistical variables will give rise to a reduction of the system's complexity. 

Except for an overall scale constant and a convenient re-scaling of variables, the information metric on the 21- 
dimensional embedded manifold A4s is given by (see Appendix A for the explicit derivation), 



dS 2 Ms 



E 



(d^lj-i + 2r 2j ^idfi 2j _ 1 dcr 2j ^ 1 



2du%_ x ) 



(6) 



where the coefficients r 2 j-i are defined as 



def 



dfi 2j dy 2j 



(7) 



The explicit expressions of such coefficients depend on the functional parametric form given to the embedding con- 
straints fi 2 j = ^2j \P"2j—ii a 2j-i)' More details are given in Appendix A. From ([7]) it follows that the coefficients r 2 j-i 
are non-zero if and only if fi 2 j depends on both (J- 2 j- X and o~ 2 j- X . Therefore, we may conclude that the emergence of 
the non-trivial off-diagonal terms in ([6]) is a consequence of the previously mentioned correlational structure arising 
from the embedding constraints. Motivated by these considerations, we will name from now on the coefficients ?*2j_i, 
macroscopic correlational coefficients. 
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B. Information Geometry of the Model 



In this Subsection, we discuss the main steps leading to computation of the asymptotic temporal behavior of 
the dynamical complexity of geodesic trajectories for the 2Z-dimensional Gaussian statistical model. For the sake 
of notational simplicity and in view of the involved analysis that we will present shortly, we replace in @ with 
9 = lr .., n u a x ,..., at) and P (x\Q^ in @ with P (X\S) so that M ^ mbcddcd ) = Ms _ {p(X\Q)}. We begin to 
study the IGAC arising from the Fisher-Rao metric defined as, 



a 2 — f 
as M s ~ 9ij 



(9) rfO'dO-' = (—d/i 2 . H ^-dfifrdak H — ^da\ ) ' w ith j = !>•-•) 



fe=i 



2/. 



(8) 



We assume positive macroscopic correlational coefficients r k £ (0, 1), Vfc = 1,..., I. The Fisher-Rao metric tensor 



(O) = gij (^j,...^; cti,..., a;) leading to the line element in © is given by, 



[<fe (9)] 



21x21 



( Mi% 



V o 





■ 

• 

M% 2 J 



\ 




(0 



with i, j = 1,..., 21. 



(9) 



where M^ 2 is the two-dimensional matrix defined as, 



M. 



(fc) 



del' 



1 rfc 

r fc 2 



with fc = 1,..., I. 



(10) 



The inverse matrix 



M. 



(fc) 



useful for computing the Christoffel connection coefficients and other quantities char- 



acterizing the information geometry of M s is given by, 



M. 



(*) 



- 1 dcf a\ 



2 -r* 
-r* 1 



with fc = 1,..., I. 



It can be shown [43| that the scalar curvature of such 2Z-dimensional manifold is given by, 

i 

K Ms (r lv .., r i )--2^(2-r2)" 1 . 



(11) 



(12) 



fc=i 



Notice that in the limit of vanishing coefficients {r k = 0}, TZm, = —I a s shown in [4J|. The computation of geodesic 
equations on the 2Z-dimensional Gaussian statistical manifold Ai s leads to the following coupled systems of nonlinear 
second order ordinary differential equations, 











d 2 Vk 
dr 2 

d 2 a k 
dr 2 



1 



r\ a k 



1 1 

2-r 2 a k 



dr 

dfJ-k 
dr 



4 1 d/j, k da k 



2r k 1 



2 — r k Ok dr dr 2 — r 2 a k 



dak 
dr 



2r k 1 d[i k da k 2r 2 - 2 1 / da k 



2 — r\ a k dr dr 



2 - r{ Ok \ dr 



(13) 



with fc = 1,..., Z. When r k — > 0, Vfc we get the ordinary Gaussian system of nonlinear and coupled ordinary differential 
equations. Integration of such coupled system of nonlinear second order ordinary differential equations in is highly non 
trivial. However, this problem can be tackled using the information geometric diagonalization procedure introduced 

in [43j]. The information metric tensor g (^ 1; ...^ ; ; ai,..., a{) = f g (0) in (jHJ is symmetric and therefore diagonalizable. 
The eigenvalues of such matrix are, 



ot± (r k ) = z , A (r-fe) 



1 + Ar%, with fc = 1,..., I. 



(14) 



The eigenvectors = O+ (r k ) and 0_ = f 6_ (r) corresponding to a+ (rfe) and a_ (r k ), respectively, are, 



©+ (r fc ) 



l+yggg J and, 9_ (r fc ) 

2r fc / 



2r fc 



with fc 



1,..., Z. 



(15) 
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The diagonalized information matrix 



g> (e (e 



\g(e)] B<M = Guxii (n,..,n)[s' (e (e 

where, in an explicit way, we obtain 



in the new basis B ncw satisfies the following relation, 



E 2lx2l (n,-. n): 



e 



^2X2 \ 

-00 

0-0 

V D® 2 J 



(k) 

with the two-dimensional diagonal matrices D 2> ! 2 defined as 

1 



D 



del' 



2 



<(^ **) \ o 3 +V^) 



with fc = 1,..., Z. 



(16) 



(17) 



(18) 



The columns of the matrix E 2 i X 2i (n.,---> n) encode the eigenvectors of [g (@)] Bold , @+ and 8_ and is given by, 



[E(r u ..., ri )} 



21x21 



where the two-dimensional matrices E 2 k ^ 2 are, 



{ E$ 2 \ 
-00 
0- 
V E® 2 ) 



(19) 



E. 



(k) 

2x2 



def 



l-y/A(r k ) 1+ V /A(r fc ) 
2r k 2r h 



(20) 



The relevance of E 2 i x2 i (ri,..., r;) (and its inverse) is in expressing the set of macrovariables (/i 1; .../i ; ; cti,..., ct/) in 
terms of the new statistical variables (fi lv .., cti,..., ct/), 



From differential geometry arguments (45|, it follows that 



= -E2ix2i (ri,..., n) 



da 



and, 



Mi 



= #2ix2i (n, — ; r fc) 



/'l 



a-; 



(21) 



(22) 



Substituting (|20)) in (1221) , we finally obtain the formal relation between the old and new set of macrovariables labelling 
the 2/-dimensional macrostates of the embedded Gaussian statistical model in presence of non-trivial off-diagonal 
terms, 

/- -vdef_,_ , - > def 1 - \/A(r fc ) . 1 + y/Ajrk) . , , 

Mfc U"*, o-fc) = Mfc + o-fe and, <r (/x fe , er fc ) = fi k H a kl with fc = 1,..., I. (23) 



2r fc 



2r k 



After having introduced the information geometric diagonalization procedure, the new line element 
ds' 2 (fa,..., fa; (Ti,..., &i) to be considered becomes, 



ds' 2 (/*!,..., fa; ax,..., ai) = ^ 



k=i L 



a- (r k ) 



rdp.1 



dal 



2 r^k ' 2 k 

[a (r k )Jl k + a\ (rk) &k] [ao (r k ) jl k + ai (r k ) a k 



(24) 



where, 



def 3 ± V A / \ def 1 V A ( r fc) / s def 1 + V A ( r k) , A , s def 1 , . 2 ,„_, 

a± (He) = o ' °o ( r fc) = o > °i ( r fc) = « an d, A (r fc ) = 1 + 4r fe . (25) 



2r fc 



2r fe 
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Notice that ds (fa,..., fa, &x,..., ai) can be rewritten as, 



ds' 2 (fa,., fa, ai,., cr/) = 



fc=i 



a- (r k ) 1 



[ffll (^fc)] <5fc fl , £o(rfc) Afe A 



As a working hypothesis, we assume that a °( rk \ ^ k ^ T \ « 1 for r > 1 and for each k = 1,.., Z. Stated otherwise, we 
assume that 



(rk)] ^(i I °o(^) fifcV 



rda? 



(26) 



lim 



Afe ( T ) 



CTfc (t) 



<C mm 

rS(0,l) 



a (r-fc) 



mm 

i- fc e(o,i) 



1- 



2.6. 



(27) 



Then, in the asymptotic long-time limit [46?], the notion of distinguishability between probability distributions on the 
diagonalized statistical manifold is quantified by the following line element, 



ds' (fa,-, fa, ai,-, crj) 



a- (r k ) 1 



[ai (rk)] °k 



12* 



r,2 



Q+ (r-fc) 1 
[ffll (nfc)f °fe 



Recall that the Christoffel connection coefficients T^- are defined as, 



(28) 



(29) 



Substituting the metric tensor components from (|28p into ([2U|) . it turns out that the only non-zero connection coeffi- 
cients are given by, 



o-fc 



fe _ a_ (r fc ) f 



i 



(30) 



where k = 1,.., Z. Therefore the set of coupled nonlinear ordinary differential equations satisfied by the geodesic 
trajectories becomes, 



d fa 2 dfa da k = Q d 2 a k | q_ (r fc ) 1 / d/2 fc 
dr 2 a dr dr ' dr 2 a+ (r^) er \ dr 



a-(rfc) 



= 0. 



(31) 



a+ ^ rfc ) — ^ ^ and the system of equations pip describing the asymptotic behavior of 
maximally probable trajectories on the diagonalized manifold becomes the standard two-dimensional Gaussian system 



Notice that in the limit of r k —> 0, 
maximally probable trajectories oi 
of nonlinear coupled ordinary differential equations studied in [36| . In order to further simplify the integration of (|31[) , 
consider the following (invertible) change of variables, 



Substituting 



r~ - \ ~ \ /2a_(rfe)_ , _ . _ \ 

(fx k , cr fe ) — > \[M k (fx k , cr fe ) = W -^"T^y^fc, ^ (Mfe: o-fe) = o-fc I • 

into (f3"Tj) . the coupled system of nonlinear differential equations to be integrated becomes, 



d 2 fa 2 dfa da\ 



d 2 a[ 



1 



dr 2 



a' dr dr ~ ' ' ^ ' I w^- 



dr 2 2cr 



dr 



1 /da', 



dr 



= 0. 



Integrating (j3"3"|) leads to the following geodesic trajectories, 

1 



fa\(r) 



"2 
"fc 



2A t 



exp(-2A/tr) + 



4A fe , <7 fc (r) 



^fc exp (-Afc-r) 
exp (-2Afcr) + 



where and Afe are real and positive constants of integration [36]. Using (|23p and (|32p . we have 



(32) 



(33) 



(34) 



Mfe (Mfe, o-fe) = 



, dcf / a+ (r fc ) , 



2a- (r fc ) 



Mfe + o fc and, a k (fa, a' k ) 



dcf 1 - y/A (r- fc ) a+ (r fe ) , 1 + ^ (rfc) , 



2rfe 



2a_ (r fc ) 



Mfe 



2r fc 



(35) 
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Notice that our working hypothesis (|27l) is satisfied since we have, 

Afc (t) 



lim 

t->oo CTfc [T) 



— • hm cx exp -A fc r — >• 0. 

2a_ (r fe ) / r^oo cr^ (r) 



Finally, in terms of the original macrovariables (fi k , Ufc), the geodesic trajectories become, 



Mfe (t; 



ot+(rk) 



0-fc(.T, Tfcj - 2rfc A/ 2a_(r fc ) 



2Afc exp(-2A fc T) + -=l 



2Afc exp(-2A fe r) + ^- 



4Al 



4A fc 



Hfc exp(-A fc r) 

32~ j 

exp(-2A fc r) + -=4- 



l + % /A(r fc ) 5 fc cxp(-A fc r) 



oxp(-2A fc r) + i ^- 



(36) 



(37) 



In our probabilistic macroscopic approach to dynamics, the geodesic trajectories in (|37p represent the maximum 
probability paths on the 2Z-dimensional embedded Gaussian statistical model. 



IV. INFORMATION GEOMETRIC COMPLEXITY 



In our information geometric approach a relevant quantity that can be useful to study the degree of complex- 
ity characterizing information-constrained dynamical models is the information geometrodynamical entropy Sm 3 ( t ) 
(IGE) 36]. In what follows, we will briefly highlight the key-points leading to the construction of such quantity. 

The elements (or points) {P (X\Q)} of a 2Z-dimensional curved statistical manifold M. s are parametrized using 

2Z-real valued variables 



1)2/ 



m s d = {p(x\q) : e = (V,..., e 21 ) g ?i tot) } 



(38) 



The set T>q°^ is the entire parameter space (at disposal) and it is a subset of 



i:2/ 



V, 



(tot) def 



)Igk = (Xgl ®Iq2...®Iq2i) C 



(39) 



fe=i 



where Igk is a subset of R and represents the entire range of allowable values for the macrovariable 9 . For instance, 
considering the statistical manifold of one-dimensional Gaussian probability distributions parametrized as usual in 
terms of = (/z, a), we obtain 



V, 



(tot) 



= 1^ ®Z a = [(— oo, +oo)(g)(0, +oo)]C 



(40) 



In the IGAC, we are interested in a probabilistic description of the evolution of a given system in terms of its 



correspondent probability distribution on Ai s which is homeomorphic to T>, 



(tot) 



Assume we are interested in the 



evolution from r^tia] to rfi na i. Within the probabilistic description, this turns out to be equivalent to study the 
shortest path (or, in terms of the ME method |27H3G |. the maximally probable path) leading to Q (rg na i) from 

© (Tinitiai)- 

Is there a way to quantify the "complexity" of such path?. We have proposed that the IGE Sm s ( t ) is a good 
complexity quantifier (25l . l26j . A suitable indicator of temporal complexity within the IGAC framework is provided 
by the information geometric entropy (IGE) Sm„ ( t ) 136| . 



S Ma (t)= tog vol \ V ^ codcsic) (r) 



(41) 



The average dynamical statistical volume vol 



(geodesic) 

e 



(r) 



geodesic) 

e 



= f lim ( — I dr' vol 



is defined as, 

(geodesic) ,f\ 



(42) 
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where the "tilde" symbol denotes the operation of temporal average. The volume vol 



(geodesic) 

e 



is given by, 



vol 



(geodesic) / / 



(r' 



del 



P{Msig) {e\...,6 n )<t®, 

iSlc V') 



(43) 



where p<M a , g) iP l >•••! ^™) i s the so-called Fisher density and is equal to the square root of the determinant g 
|det {g^u)\ of the metric tensor g^ (0), 



P{M„ a) 



del' 



A ((^,-,0)1 



(44) 



The integration space 2?(s codcslc ) (7"') in (|43|) is defined as follows, 



(geodesic) 



(t') d = {e = ((9 1 ,..., 9 n ) : 9 k (O) < fe < 6> fe (r')} 



(45) 



where k = 1,.., 71 and 9 =9 (s) with < s < r' such that, 



ds 2 



c?s ds 



(46) 



The integration space 

pfeeodesie) ^ ^ 

is a 2^-dimensional subspace of the whole (permitted) parameter space 
2?g ot) . The elements of<° dcsk) (r') are the 2Z-dimensional macrovariables {0} whose components 9 are bounded 
by specified limits of integration 9 k (0) and 9 k (r') with k = 1,.., 21. The limits of integration are obtained via 
integration of the 2/-dimensional set of coupled nonlinear second order ordinary differential equations characterizing 
the geodesic equations. Formally, the IGE Sm s (t) is defined in terms of an averaged parametric (r is the parameter) 
21 + 1-fold integral over the multidimensional geodesic paths connecting (0) to (r). In our information geometric 



approach, the information geometric complexity vol 



(geodesic) 

e 



(r) 



represents a statistical measure of complexity 



del 



of the macroscopic path = (r) on Ais connecting the initial and final macrostates 0; and 0/, respectively. 
The path (r) is obtained via integration of the geodesic equation on Ms generated by the universal ME updating 
method. At a discrete level, the path (r) can be described in terms of an infinite continuos sequence of intermediate 
macroscopic states, (r) = [©»,..., 0|_X) ®k> ©fe+iv) ©/] with Qj = (tj), determined via the logarithmic relative 
entropy maximization procedure subjected to well-specified normalization and information constraints. The nature of 
such constraints defines the (correlational) structure of the underlying probability distribution on the particular curved 
statistical manifold Ms- In other words, the correlational structure that may emerge into our information-geometric 
statistical models has its origin in the valuable information about the microscopic degrees of freedom of the actual 
physical systems. It emerges in the ME maximization procedure via integration of the geodesic equations defined 
on Ms and it is finally quantified in terms of the intuitive notion of volume growth via the information geometric 
complexity vol 2?te codcslc ) ( r ) 0T ^ m entropic terms by the IGE Sm s ( t )- The information geometric complexity is 

then interpreted as the volume of the statistical macrospace explored in the asymptotic limit by the system in its 
complex evolution from 0, to 0^-. Otherwise, upon a suitable normalization procedure that makes the information 
geometric complexity an adimensional quantity, it represents the number of accessible macrostates (with coordinates 
living in the accessed parameter space 2?te eodcslc ) (7-)) explored by the system in its evolution from 0^ to 0/. 
For the model defined in ((SJ, Vm s ( r ) becomes, 



Vm s (t) = vol 



V, 



(geodesic) 

e 



(r) 



d ^ f i Idr' 



(47) 



e«(o) 



where g = det [gij (0)] 2 / X 2Z ^ s ^ ne determinant of the block-diagonal matrix in ([5]), 



■(r 1 ,..,r ; ) = n det ( M 2x ) 2 )=n 



fc=l 



fc=l 



(48) 
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The geodesic paths (r) = (n 1 (t; rx) , cti (t; ri) ;..; [i l (r; n) , 07 (r; n)) are given in (l37l) . Substituting them into 
63), we get, 



V Me (r; Ax,.., X k ) = J] ^ ^ 



fc=i 



S fc exp (-AfcT 7 ) - 4A fc ^/^t(^ y exp (-2A fe r' 



— ^ n fc exp (-A fe r') - 4A fc w 2a + _y fc ) v 2rfc exp (-2A fc r') 



dr 1 



For the sake of simplicity, let us introduce the following substitutions, 



. dcf „ D dcf . -. 



a+(rfe) dcf 1 + V A ( r fc)~ n dcf . 

, <^k = 7. ^k, L>k — — 4Afc 1 



, 2a_(r fc )' " 2r fc 
Then, the integral defining Vm s ( t ; Ai,..,Afc) becomes 

/ 



/ a+ (r fc ) 1 - y/Afa) 
2a_ (r fc ) 2r fc 



(49) 
(50) 



V.m s (t; Ai,..,A fe ) = ]J[ 



fc=X 



A 



- XkT ' ' B k e- 2XkT ' 



C k e~ XkT ' +D k e- 2X " T 



rdr' 



dr' 



(51) 



Upon integration, we get 
he 



T A- "~ XkT ' 



C fe e~ A ^'+^ e -2A fe r' A fc VCfe ^fe 

and substituting (|52|) into f|5 1 [) . we obtain 

i 



D k + C fe e A '= T 



^ fc r r ~°° - ( ^ k ^ k \\n^ k t (52) 
Cfc Afc \C'k Dk J Dk Ck 



V Me (r; A 1 ,..,A fe ) = n{y / 2 _ 

fc=i 



C k Xk \C k Dk) n Dk t 



(53) 



Introducing again the original parameters in (|50[) . we finally get 
(r; A lv .,A fc ) = JJ J i 



V 



Ms 



1 1 1 + v%) 



2r fe V2" 



2r fc ^ 



1 + v^Mj A fe (l - vaw) A fe 

where the strictly positive function E (rfc, Afc, a±) is given by, 

l+VA(r fc ) 



InE (rjfe, Afe, a±) 



(54) 



E (rfc, Afe, ay.) 



dcf 



4\, / 1 ~\/ A ('-fc) 

^M/ 2a„(r fc ) 2r fc 



Finally, inserting (|54|) into (|4T|) . the IGE <Sa/i s (r) becomes, 



Sm. (t; Ai,..,Afe) = ^ log ^ - 

fe=i 



2r fev ^ 



2r fev /2^! 



> 0, Vr fc e (0, 1) . 



2r fcv ^7f 



(55) 



In E (r fc , Afc, a±) 



1 + VAW) Afc (l - VA(^)j Afc _ 
With a suitable change of notation, equation (151))) can be rewritten in a more elegant way as follows, 

A 2 (rfc, Afc) 



(56) 



Sm s (r; Ai,..,A fc ) T ~°° ^log 



fc=i 



Ai (rfc) + 



where, 



Al (r fe ) d ^ A 2 (r fc , A) + ^) (2 - rjj) ln£ (r,, A <e ,a ± ) > ^ ^ dcf WT+4^ 



rfc 



Ai 



(57) 



(58) 



As stated above, E (r k , Afc, a±) is a strictly positive function of its arguments. It appears that the introduction of 
embedding constraints between the macrovariables of the larger Gaussian statistical model leads to the emergence of an 
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asymptotic information geometric compression of the explored statistical macrostates on the embedded configuration 
manifold M. s in its evolution between the initial and final macrostates. This result, thought for a special (not general) 
case, leads to interesting conclusions. The presence of constraints between macroscopic pieces of relevant information 
on the microscopic degrees of freedom of a complex system allows for an information geometric probabilistic description 
whose complexity, measured in terms of Sm 3 ( or Vm„)> decays in a power law way. Asymptotically, the complexity 
reaches a saturation value characterized solely by the strength of such macroscopic correlational coefficients. The 
relevance of such results becomes even more clear if compared to what happens for the 4i-dimensional (uncorrelated 
and larger) Gaussian statistical model [36j . We will mention this comparison in our final remarks. For = r s , Vfc, 
s = 1,..., I, the information geometric entropy Sm 3 ( r ! I, A, r ) becomes, 



S Ms (t; I, A, r) 



log 



Ai (r) + 



A 2 (r, A) 



(59) 



Therefore, the information geometric complexity presents a power law decay where the power is related to cardinality 
I of the microscopic degrees of freedom characterized by correlated pieces of macroscopic information and it reaches 
a saturation value quantified by the set of coefficients {r k }. 



V. JACOBI-LEVI-CIVITA EQUATION AND JACOBI FIELDS 

The Jacobi-Levi-Civita ( JLC) equation of geodesic deviation is a complicated second-order system of linear ordinary 
differential equations. It describes the geodesic spread on curved manifolds of a pair of nearby freely falling particles 

travelling on trajectories P (r) and /p (r) = P (r) + 5@ p (t) . The JLC equation is given by [13] , 

l^+^^i^O, (60) 

with k = 1,.., 21 and where the covariant derivatives V ® v ^ along the curve P (t) are defined as, 

V9»(t) d _e f d9"(r) d&? 

Vt dr vp dr K ' 

The Jacobi vector field components J k are given by, 

J fc ^^e^ f ( ae ^ Afc) ) *A fe , (62) 
and TZafs^s is the Riemann curvature tensor defined as [47| . 

T^lup = ®vTp p — d p T™ v + Tp v T^ p — r^ p r^. (63) 

In the 2/-dimensional case J — {J k } k _ x 2 z re P resen ts how geodesies, in a 1-parameter family of geodesies, are 
separating. The covariant derivative in (|60l) is defined as (48j . 

2? 2 J p _ d 2 J" „ d.J a dQf 3 „ a d 2 Q? M d& v d<Q p a ^ a d<d a dQ 13 p 

~ + 2r ^^7^T + T ^ J ~d^ + T ^ u ~dr~~dr~ + r ^ rp -^T^T J ■ (64) 

Equation (|60[) forms a system of 21 coupled ordinary differential equations linear in the components of the deviation 
vector field (|62|) but nonlinear in derivatives of the metric tensor gij (0). It describes the linearized geodesic flow: 
the linearization ignores the relative velocity of the geodesies. When the geodesies are neighboring but their relative 
velocity is arbitrary, the corresponding geodesic deviation equation is the so-called generalized Jacobi equation [49| . 
The nonlinearity is due to the existence of velocity-dependent terms in the system. Neighboring geodesies accelerate 
relative to each other with a rate directly measured by the curvature tensor TZ a p^s- 

Multiplying both sides of (|60[) by (0) and using the standard symmetry properties of the Riemann curvature 
tensor, the geodesic deviation equation becomes, 

V 2 J l m dQ k rm dQ l n 
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After some algebra, it follows that the only non-zero Riemann tensor components are given by, 



(V \ k - a ~ ^ 1 ft- - 1 91 



(66) 



In the model considered, the Jacobi held J» has 2Z-components J* 1 that can be grouped into Z-pairs as 

I J fl=l,.., 21 

follows, 



1-pair 



r2fc-l J2k 



fc-pair 



( 1 % 7 X 



Z-pair 



, (67) 



where, 



m-i def / dji k (r; A fc ) 



x . J2k def / ggfc (r; Afc) 

6xk ' JM ° - ( d\ k 



5 At 



(68) 



with fc = 1,.., /. Similarly, the 2Z-equations of Jacobi-Levi-Civita can be grouped into Z-pairs of differential equations 
with identical structure. Substituting (|3U|) . (j6"6")l and into ([63)1 , after some tedious algebra, the hrst pair of JLC 
equations become, 



_ d 2 J 2k ^ , x ,k db k dj 2k ^ k dji k dJ 2k 



d T 2 



dr dr 



?2fc-l 



(ria) 



dr 2 



i\t/ r i\' . (-R1212) 



dr 



j2k 



[ 12 > ~d^ 



d*> (ri 2 ) 



/ r i \*/Vii \ fc (-R1212) \ dji k da k 

[Ll2) [Ll2 > > ^ 1 "d7^7 



(5n)" 



(69) 



and, 



2 t2/c 



d 2 J 



dT 2 



. 2(r?i) ^^ +2(ri2 )' 



dr dr 



2fe-l 



(r 11 )' s 



fc d 2 fl k 



dr 2 



fc. (HO* + (r?0* (r} 2 ) fc + (H 2 ) fc (r?0* - (iW) » d ~° k 



(922) 



dr dr 



(r| 2 ) 



fe d 2 a-fc 
"dr 2 " 



fc* (ry fc + (ry fc (ry fc ) 



(322) 



.(70) 



More explicitly, equations (|69|) and ([70]) may be written as, 



= 



d 2 J 2k ~ l 2 da k dJ 2k ~ x 2 djl k dJ 



2k 



-J 



dr 2 b k dr dr a k dr dr 
2k \ 1 rf2 Afc , 3 dp, k da k 



+ J 



2fc-l 



1 did 1 I d<7, x ' 



<5fc dr 2 ai 



dr 



<j k dr 2 a\ dr dr 



(71) 



and, 



d 2 J 2k n 



a_ (r fc ) 1 d/2 fe dJ 2fc - 1 2 d5- fe dJ 2fc 
a+ (rfe) CTfc dr dr a k dr dr 



j2k-\ 



a-(r k ) 1 d 2 fl k a_ (r k ) 1 dfi k da k 



a+ (r k ) a k dr 2 a+ (r k ) a\ dr dr 



j2k 



1 d 2 a k 2 f da 



Ok dT 2 a 2 \ dr 



- 2 



a_ (r k ) 1 / dji k 



a + i r k) o\ V dr 



(72) 
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From ([32]) and (p4|) . we notice that the asymptotic expansion of a k and the first and second derivative of jl k and a k 
are given by, 



(r) 
dr 2 



8At, . . d/ij, 
exp (-A fe r) , -p 



1 a+ (r k ) 64\j da k 8A£ 

^-^2" exp ~2A fe T , -7— « — exp (-A fe r) , 

2a_ (r fc ) ^ 



a+(r fc ) 128A* . , . d 2 a fc _ 8A^ 



2a_(r fc ) £ 2 



*L e xp(-2A fe T), 



dr 2 



^exp(-A fc r) . 



(73) 



Substituting ([73)) into ([69)) and ([70l) , keeping only the leading terms into the asymptotic expansion, the JLC equations 
to integrate become, 



d 2j2k-l d j2k-l 

+ 2\ k - 



dr 2 



and, 



dr 



2 j2k 



1 a+ (r k ) IQK 
2a_ (r k ) E k 



exp (-A fe r) 



dJ 



2k 



d 2 J 



1 2a- (r k )8\ 2 k 



dr 2 



a + (r k ) S fe 



7=^ exp (-A fc r) 



dJ 



dr 



2k-l 



1 a+ (rk) 8\j 
2a- (r fe ) S fe 



exp(-A fe r) J 2fc = 0, (74) 



dr 



2A, 



dJ 



2 k 



dr 



A 2 .J 2fe = 0. 



As a working hypothesis, we assume that [36J, 

dJ 2k ~ x 



lim 

r — »oc 



exp (-A fe r) ■ 



dr 



0, lim 



exp(-A fc r) 



dJ 



2 k 



dr 



0, lim 

T— t-OO 



exp (-A fe r) J 



21. 



0. 



(75) 



(76) 



In order to prove that our assumptions in (|76[) are correct, we will check a posteriori their consistency. The geodesic 
deviation equations in ([74]) and ([751) finally become, 



2 ?2fc-l 



d 2 J 



dJ 



2fc-l 



2 j2k 



d 2 .J 



dJ 



2k 



2\ k — = 0, -tV + 2A fc - 7 — + \ 2 J 2k = 0. 



dr 2 g?t dr 2 * dr 

Integration of ([77l) leads to the following asymptotic expressions for J^ 1 (t) and J 2 ^ (t), 

J^" 1 ( T ) = C o fe) + Cf > exp (-2A fc r) , J 2k s (r) = cf 5 exp (-A fc r) + Cf V exp (-2A fe r) 



(77) 



(78) 



where Cw" 1 are reed constant of integration with w = 1,.., 3. Notice that conditions (|76p are satisfied and therefore 
our assumption are compatible with the solutions obtained. 

Consider the Jacobi vector field components w M f defined in ([67| and its magnitude Jm s , 

I ) ii— i,.., i 



j2 dc 



The magnitude Jm s is called the Jacobi field intensity. Using ([28]) . Equation (|79[) becomes, 

i 



a- (r k ) 1 

it 1 [fll ( r fc)f ^fe 



2 q+ (^fc) i 

[ai (r fe )] 2 CTfc 



(79) 



(80) 



Substituting ([75)) and ([75)) into (1501) . keeping only the leading terms in the asymptotic expansion of J 2 ^ (r), we 
obtain 



3U. (r) ™ £ 



fc=i 



_o 

8A? 



exp (2A fe r) 



(81) 



Let us rewrite the quantity J 2 ^ (r) in terms of "attenuation factors" A(r k ) and elementary quadratic Jacobi vector 
field components j 2 M (r; Afe) given by, 



drf a_ (r fe ) 



2r fc (3 - yiTlr 2 



(82) 
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and, 



3m, ( r ; X k) 



( n( k )~ 
def / 'k 



exp (2A fc r) 



(83) 



respectively. Notice that A k (r k ) is a bounded function of the correlation coefficient r k £ (0, 1) and its maximum is 

(84) 



reached for r k = y/2 — \/2 ~ 0.77, 

Amax (r-fc) = f max A k (r k ) =A k (r k ) = 3- 2^2 ~ 0.17. 



r k e(0, 1) 

Therefore, substituting (|82|l and ([S3]) into (f5T|) and considering the boundedness of the attenuation factors in 
the square of the Jacobi field intensity J 2 M (r) may be written as, 



J 2 M 3 ( r ; A lr-: ^l) T «°° A {rk)3M 3 ( T 5 A fe) • 



fc=i 



Let us now consider the asymptotic behavior of Jm 3 j 



7 2 def 
J A4 S — 



fc=l 



SB' 



J 



(2fc-l) 



2^fc t-(2A-X) 7 (2/c) _2_ 
9 J M a J M S + „2 



7(2fc) 

y M s 



(85) 



(86) 



From ([23)1 and (|62p . it follows that the Jacobi field components {j^, 1 \ Jjw ^ with fc = 1,.., Z are related to the 

(~(2k—l) ~(2k)\ 
^m, ' JjA s ) m tne following wa Y: 



7 (2fe-l) def ?(2fe-l) T(2fe) 



(87) 



and, 



T (2k) def 1 ~ x/A(r fc ) ~( 2 fc_i) 1 + \/A (r fc ) ~( 2 fc) 
J ^ " 2r fc ^ + 2r, ^ ' 

Considering (|78|) and substituting (|87| and (|88|) into (j86|) . the asymptotic behavior of is given by, 



J 



M, 



£ 

fc=i 



2r fc 



2r fc 











i 


" f(2fc-l)" 







(89) 



Substituting (f3"T|) and (|75|) into and keeping only the leading terms into the asymptotic expansion, we obtain, 



(40 



embedded 



(t; Ai,.., A;) fa A k {r k )j 2 M 3 ( T ; A fe) > 



k=l 



where 2m ( r ; Afe) = 2m (r; A&) is defined in (JSHJ and the new attenuation function is given by, 

def 

4r-2 [i + 2 r fe a (r fe )+a^ (r k )] 



A k (r k ) 



[1 + oi (r fc )]' 



(90) 



(91) 



where ao (r^) and ai (r^.) are defined in (I25p . Notice that _4fc (r^) is a bounded function of the correlation coefficient 
r k £ (0, 1) and its maximum is reached for f k ~ 0.65, 



def 



-Afcmax (r k ) = max A k (r k ) = A k (r k ) ~ 0.15. 

i-fe£(0, 1) 



(92) 



In [36j,|44|, it was shown that in absence of constraints, the model here considered leads to an asymptotic behavior of 
the Jacobi fields given by, 



(93) 



k=l 
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Therefore, from equations (|93l) and (1901) . we obtain, 



< 



( J aU embedded A l'"> A ') 



V 



— k ' ov kn < 0.4 < 1, 94 

[l + «i(r S )] 2 



We conclude that the appearance of embedding constraints among the Gaussian statistical macrovariables on the larger 
4Z-dimensional curved manifold leads to an attenuation of the asymptotic exponential divergence of the Jacobi field 
intensity on the embedded 2^-dimcnsional manifold. This is a quantitative indication that the information geometric 
complexity of a system decreases in the presence of emerging correlational structures. 

VI. FINAL REMARKS 

In this article, we characterized the complexity of geodesic paths on a curved statistical manifold M. s through the 
asymptotic computation of the information geometric complexity Vm s and the Jacobi vector field intensity Jm 3 - We 
considered a manifold a 2Z-dimensional Gaussian model A4 S reproduced by an appropriate embedding in a larger 
4£-dimensional Gaussian manifold and endowed with a Fisher-Rao information metric v (©) with non-trivial off 
diagonal terms. Such terms in the information metric on the embedded manifold emerged due to the presence of 
a correlational structure (embedding constraints) among the statistical variables on the larger manifold and were 
characterized by macroscopic correlational coefficients r^. First, we observed a power law decay of the information 
geometric complexity at a rate determined by the coefficients and concluded that the non-trivial off diagonal 
terms lead to the emergence of an asymptotic information geometric compression of the explored macrostates on 
Ai s . Finally, we observed that the presence of such embedding constraints lead to an attenuation of the asymptotic 
exponential divergence of the Jacobi vector field intensity. 

The relevance of such results becomes evident when compared to what happens in the larger Gaussian statistical 
model (absence of constraints) [H, Hi]. In such case, the information geometric entropy of the 4/-dimensional (larger 
and uncorrelated) Gaussian model increases linearly in time and its complexity diverges exponentially at a rate 
determined by Afe, the Lyapunov exponents of the statistical trajectories of the system [25j and the Jacobi field 
intensity diverges exponentially without any attenuation factor. It seems that our measure of complexity not only 
can quantify the degree of chaoticity of a physic al sy stem, but it also adequately captures the correlational structure 
(relationship between system's components, [al9l-[lllp in its behavior. In the model studied, the emergence of structure 
appears in terms of non-trivial off diagonal elements in the Fisher-Rao information metric. Such structure leads to 



£)(gcodcsic) ^ 



and, thus, to a reduction of the complexity of the path 



the information geometric compression of vol 

leading to fina i from ©initial- 

Information Geometry and Maximum (relative) Entropy methods hold great promise for solving computational 
problems of interest in classical and quantum physics in terms of their probabilistic description on curved statistical 
manifolds. Our theoretical formalism allows us to tackle physics problems through statistical inference and informa- 
tion geometric techniques, that is Riemannian geometric techniques applied to probability theory. The macroscopic 
behavior of an arbitrary complex system is a consequence of the underlying statistical structure of the microscopic 
degrees of freedom of the system being considered. 

As a side remark, we point out two more facts: 1) Probabilistic concepts are naturally incorporated into the 
fundamentally statistical quantum theory. Furthermore, describing and understanding the complexity of quantum 
motion is still an open problem since our present knowledge on the relations among complexity, chaoticity and 
dynamical stochasticity are not satisfactory at all 2) Riemannian geometric tools are currently being used to 
characterize the quantum gate complexity in quantum computing [5(il - l5ii ] . In (Hoj , the problem of finding quantum 
circuits was recasted as a geometric problem. It was shown that finding optimal quantum circuits is essentially 
equivalent to finding the shortest path (geodesic) between two points in a certain curved geometry. In light of these 
two considerations and in view of the results obtained thus far, we are confident the work presented here constitutes 
a further important step towards the characterization of the dynamical complexity of microscopically correlated 
multidimensional Gaussian statistical models, and other models of relevance in more realistic physical systems. We 
hope to extend this approach in the field of Quantum Information to better understand the connection between 
quantum entanglement and quantum complexity [5, 

mil, in in. 
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Appendix A: Derivation of the line element 

We derive Equation (|6]). For the sake of clarity, we consider a two-dimensional Gaussian probability distribution 
p (xi, x 2 1 Hit cr i) obtained from the Gaussian distribution p (xi, x 2 \ Hi, M2, a ii <J 2) setting <Ji = a 2 an d defined as, 



p{xi, x 2 \hi, H2-, 



1 



2no{ 



exp 



1 

'2^1 



(xi - Hif + O2 - H2? 



(Al) 



The Fisher- Rao information metric in the three-dimensional statistical manifold H21 a i) ^ s given by, 



ds 2 = — (dnl + dfil + idal) 



(A2) 



Consider the two-dimensional submanifold embedded as a slice in the three-dimensional space defined by the following 
embedding constraint, 



M2 = Ma (Mi, ■ 



(A3) 



The Gaussian distributions p(xi, X2\Hi> M2 (Mi, °"i) 1 a i) = p(xi, x 2 \hi, &i) belonging to this submanifold are such 
that (X2) = H2 an d = Mi are n °t independent, they are related in a peculiar way. From (|A3I) . we obtain 



, &H2 , ^M2 , 

dfi 2 = ^ — dHi + ^ — 



9 Hi 



dai 



(A4) 



that is, 



dnl 



dH2 
9hi 



2 dHl + (^ 2 
\oai 



daj + 2^^dHidai 
OHi 001 



(A5) 



Substituting (IA5|) in (IA2|) . the information metric becomes, 



ds 2 = — [A^^dnl + 2A l ^ iai dfi 1 da 1 + 2A aiai daf\ 



(A6) 



where the coefficients A^ 1 ^ 1 , A^ iCri and A airTl are given by, 

2 



.4 



clef 



1 



dH2 
8hi 



A 



def dH 2 dH2 a , def , 1 / 0h 2 

OHi ocr 1 2 \o<Ji 



(A7) 



Re-scaling the variables in such a way that Mi = f Ajj. ^ Hi an d ^1 = f A3 icri <7i and assuming that the coefficients 
^Mi^i an d A-(j \o \ a re constants, the line element in (|A6[) becomes 



ds 2 = A ai(7l — [dfj, 2 + 2ridfj,idai + 2da 2 ] , 
&i 



(A8) 



with, 
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From (|A7p . it follows that the embedding dehning the two-dimensional submanifold in the larger three-dimensional 
manifold of Gaussians parametrized by (p 1 , /i 2 , <J\) is given by, 



=Mi 



2 + 1 



(A10) 



with, 



Finally, using (|A10I) and (|A9|) . the explicit expression for ri becomes 



(1) dof 



(All) 



ri 



2+2- 



(A12) 



From (IA12|) . it is transparent that the explicit expression for n depends on the functional parametric form of the 
embedding constraint in (jAlOl) . 
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